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I. INTRODUCTION 

The field of multiple scattering of classical waves (elec- 
tromagnetic, acoustic and elastic waves, etc.) in dis- 
ordered media has revived in the eighties when 
the far-reaching analogies between the diffuse transport 
of waves and electrons in mesoscopic systems have been 
realized (see Ref. for a comprehensive review of the 
latter issue) . Using classical waves to study such phenom- 
ena as weak and strong localization, mesoscopic correla- 
tions and universal conductance fluctuations (see Refs. 
[sj-flfij for reviews) appears to be advantageous in many 
aspects: no need for low temperatures and small samples, 
better control of the experimental apparatus, possibility 
of more sensitive measurements, etc. In addition, ex- 
periments with classical waves can be readily performed 
in the linear regime, excluding interaction between scat- 
tered waves and hence simplifying the interpretation of 
experimental data, whereas the electron-electron inter- 
action is always present in the realm of mesoscopic elec- 
tronics and cannot simply be 'turned off', introducing 
significant difficulties in the theoretical model [TT| . 

The 'interaction' of classical waves, analogous in some 
sense to the electron-electron interaction, can come about 
if the waves propagate in a nonlinear medium. Such 
an interaction is a matter of scientific enquiry in the 
fields of nonlinear optics 0] > nonlinear acoustics [l3[ . 
etc. Nowadays, the latter are rapidly developing scien- 
tific disciplines on their own right, having considerable 
fundamental importance and numerous practical appli- 
cations. Various nonlinear phenomena (self-action: self- 
phase modulation and self- focusing of pulses and beams, 
harmonic generation, shock wave formation, etc.) occur 
in homogeneous nonlinear media depending on the type 
and strength of the nonlinearity. The effect of weak disor- 
der on the propagation of nonlinear waves has also been 
studied both theoretically (treating it as a weak pertur- 
bation) and experimentally (e.g., for laser beam prop- 
agation through atmospheric turbulence) |l4j . Unfortu- 
nately, only few experiments 0] have been performed up 
to now on multiple, diffuse scattering of classical waves in 
nonlinear media. Theoretical analysis of the problem is 
more advanced (see the bibliography of Ref. Hal), but is 
still insufficient to stimulate further experimental effort. 
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Meanwhile, diffuse waves seem to be a good candidate for 
a detailed study of the combined effect of disorder and 
nonlinearity (or interaction) on wave (or quantum par- 
ticle) propagation, as compared to interacting electrons 
in disordered mesoscopic samples: the strength and the 
type of nonlinear wave 'interaction' can be readily con- 
trolled and the nonlinear term in the wave equation is 
often of simple algebraic form. Despite these important 
simplifications, diffusion of classical waves in nonlinear 
media remains an involved problem and is far from being 
solved. The purpose of the present paper is to review our 
recent theoretical results 0-0] which are particularly 
susceptible to stimulate the reader's interest, further the- 
oretical analysis, and, hopefully, new experiments. 

For the sake of concreteness, we restrict our consider- 
ation to the problem of self-action of a scalar monochro- 
matic wave (frequency uj) in a nonlinear disordered 
medium. This is described by the following wave equa- 
tion: 

{ V ' - ^ t 1 + fe ( r ) + A£ ni>, t)]| E(r, t) 

= J(r,t), (1) 

where E(r,t) is the wave amplitude, c is the speed of the 
wave in the average medium, J(r,t) — Jo(r) exp(— iujt) 
is a monochromatic source term, fc(r) is the fractional 
fluctuation of the linear dielectric constant, and the non- 
linear part of the dielectric constant A£nl(i", t) depends 
on the intensity I(v,t) = \E(r,t)\ 2 of the wave. Eq. 
describes, e.g., propagation of light in a medium with 
intensity-dependent refractive index and (as the reader 
might already have noted) we adopt the 'optical' termi- 
nology from here on. Furthermore, we consider E in Eq. 

to be a complex quantity, thus neglecting the (possi- 
ble) generation of the third harmonics, and restrict our- 
selves to a weak nonlinearity (see below for the definition 
of 'weakness'). 

There are two fundamental questions to consider con- 
cerning the propagation of waves in a disordered medium 
with a weak nonlinearity First, one can ask about the ef- 
fect of nonlinearity on the phenomena known for waves in 
linear disordered media. We partially answer this ques- 
tion in Sec. IIVI where we summarize the results of cal- 
culation of the angular correlation functions of scattered 
waves and the analysis of the coherent backscattering 
cone in a nonlinear disordered medium. The second ques- 
tion is more challenging: can the weak nonlinearity give 
rise to new physical phenomena which are not present in 
the linear medium? The answer to the latter question is 
'yes' as we demonstrate in Sec.lVlwhere the instability of 
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diffuse, multiple-scattered waves in a disordered medium 
with a weak nonlinearity is considered. In order to fa- 
cilitate the presentation of our main results, we start by 
a very brief review of principal results known for waves 
in linear disordered (Sec.^) and homogeneous nonlinear 
(Sec. media. 



II. WAVES IN LINEAR DISORDERED MEDIA 

Waves propagating in a linear disordered medium are 
described by Eq. @J with Ae^i,(r,t) — 0. Despite its 
linearity, this problem is a rather involved one and it 
attracts a lot of attention (see, e.g., the articles in the 
present book and Refs. and [13 )• I n what follows, we 
are interested uniquely in the diffusion regime of wave 
propagation that is realized for ko£ » 1, where k =ui/c 
and £ is the mean free path due to disorder. To sim- 
plify still further the consideration, we assume that the 
correlation length of &(r) is much shorter than the wave- 
length A and adopt the model of the Gaussian white-noise 
disorder: (fe(r)fe(ri)) = ^i\j{k\t)8{v — ri), so that the 
scattering and the transport mean free paths coincide. 
Besides, we consider a disordered medium without ab- 
sorption (i.e. both Ss and Aenl are real), except if the 
opposite is specified explicitly. In the regions of space 
located at distances greater than £ from the sources of 
waves and the boundaries of the sample, the average in- 
tensity (I(r)) then obeys a diffusion equation |20j : 



0_ 

Of 



(/(r,t))=W 2 (I(r,t))+S(v,t), 



(2) 



where D — c£/3 is the diffusion constant and S(r,t) is 
the source term. In the following, we assume a time- 
independent source term in Eq. J5J: S(r, t) = S(r), yield- 
ing 7(r,t)=7(r). 

The average intensity (I(r)) given by Eq. is not suf- 
ficient to describe the intensity I(r) of the wave, since the 
latter exhibits large fluctuations: SI(r) = J(r) — (I(r)) ~ 
(/(r)) 0,|23. The spatial correlation function of SI(r,t) 
contains a short-range (Ar < £) but strong (~ 1) contri- 
bution |2|: 



(5/(r)«J(n)) = </(r)> (J(n)) 
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where Ar = r — n, Iq is a typical value of (J(r)), G(r, ri) 
is the Green's function of Eq. © with (d/dt) = 0, and 
the last expression in Eq. is obtained assuming that 
G(r, ri) is equal to its value in the infinite medium: 



G(r, n) = G (r, ri) = (AnDAr)^ 1 . The short-range 
contribution to the intensity correlation function @ de- 
scribes the speckle structure of the spatial intensity dis- 
tribution in disordered media, while the long-range con- 
tribution testifies that different speckle spots are cor- 
related (weakly, since k$£ 3> 1). The correlation func- 
tions given by Eqs. J3J and J3J are commonly referred to 
as and C^ 2 \ respectively. Other contributions to the 
intensity correlation function can exist: is a factor 

(k £)- 2 < 1 smaller than C< 2 ) and G<°) ~ (hoi)- 1 
is non-vanishing only if the source of waves has the size 
smaller or of the order of the wavelength A |2l| . We do 
not consider the two latter contributions here. 

Short- and long-range correlation functions of the field 
and intensity fluctuations can be also defined in the an- 
gular domain. For our purposes, it will be sufficient to 
consider only the short-range correlations C^- 1 ' . If a plane 
wave is incident on a disordered slab of thickness L ^> £ 
(slab surfaces being perpendicular to the z-axis) with 
a wave vector k a = {q Q ,fc QZ }, where q a = {k ax ,k ay }, 
the field (intensity) of the wave transmitted through the 
slab with a wave vector k(, is E(k ai k ) [7(k a ,k(,) = 
\E(k a , kb)| 2 ]. If the above experiment is repeated for 
different incoming and outgoing wave vectors and k^, 
respectively, the normalized correlation function of scat- 
tered fields 

G(k a , k^; k a , k b ) = C a b a 'b' 

= (E(k a ,k b )E*(k' a ,K)) , 5) 

(i(k a ,k b )) 1/2 (i(K,K)) 1/2 

appears to have a rather simple form : 

Caba'b' = S A q a ,Aq b F T (Aq a L), (6) 

where Ft(x) = x/sinha;, Aq a = q a — (and similarly 
for Aq;,), and Aq a £ <C 1 is assumed. 

The field correlation in reflection, when both incident 
and scattered waves are on the same side of the slab, is 

0ua 

Caba'b' - *Aq„,A qi [-F R (Aq a ) + F R (q a + q 6 + Aq a )] , (7) 
where 



F R (q) ~ (1 + 2C)- 2(1 + C)V + 



(8) 



£ ~ 1 is the extrapolation factor (£ = 2/3 in the diffu- 
sion approximation), <;f « 1, and the semi-infinite disor- 
dered medium is considered. The short-range correlation 
function of intensity fluctuations, G' 1 ', in transmission 
(reflection) geometry is obtained by taking the square of 
the absolute value of Eq. JBJ [Eq. Q]. 

Finally, a well-studied phenomenon in the linear disor- 
dered medium is the coherent backscattering , consist- 
ing in a two-fold (with respect to the incoherent, diffuse 
background) enhancement of the average scattered inten- 
sity in the direction of exact backscattering (k^ = — k a ). 
For normal incidence (q a = 0), the angular line shape of 
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the coherent backscattering cone can be approximately 
expressed using the function F R defined in Eq. |JSj: 

(±M „ F R (l/L a )+F R [(g 2 + l/Ll)^} 

(line) ~ F R (1/L a ) ' 1 1 

where q = and we have taken into account the absorp- 
tion of light in the medium by introducing the macro- 
scopic absorption length L a £. 

III. WAVES IN HOMOGENEOUS NONLINEAR 
MEDIA 

If Se(r, t) — 0, Eq. describes the propagation of 
waves in a homogeneous nonlinear medium (no scatter- 
ing) with an intensity-dependent dielectric constant |12| . 
Below we consider the case of Kerr nonlinearity. If the 
nonlinear response of the medium is instantaneous and 
local, the dependence of A£nl(i", i) on J(r, t) is rather 
simple: A£nl(i", t) = £2-f(r, t), where £2 is a nonlinear 
coefficient. In general, however, the nonlinear response is 
not instantaneous and can be modeled by a phenomeno- 
logical equation of Debye type for Aenl: 

d 

r NL — Ae NL (r,i) = -Ae NL (r, t) + e 2 I{r, t), (10) 
ot 

where t^l is the response time of the nonlinearity. 
Moreover, the non- locality of the nonlinear response of 
the medium can be taken into account by considering 
A£nl(f,£) obeying an equation of diffusion type: 

d 

r NL — Ae NL (r,i) = (a NL ) 2 V 2 [Ae NL (r,i)] 

- Ae Nh (r,t)+e 2 I{r,t), (11) 

where onl is some characteristic length, describing the 
degree of non-locality of the nonlinear response. 

The physics of nonlinear optical phenomena in Kerr 
media is rather rich [T^ and it is not our purpose to re- 
view it here. In the framework of our study, we would 
like, however, to call the attention of the reader to the 
following two simple but illustrative examples. First, if 
a nonlinear Kerr medium, illuminated by a monochro- 
matic plane wave of intensity Iq, is put in an optical 
resonator, the feedback provided by the resonator leads 
to the instability of the steady-state solution E(r, t) = 
Eo(r) exp(— iu>t) and the intensity I — \E\ of the wave 
develops a time dependence 22]. This phenomenon oc- 
curs only at Iq exceeding some threshold and the dy- 
namics of /(r, t) becomes progressively more complex as 
Jo increases. The second example concerns two counter- 
propagating light waves in a Kerr medium. Again, the 
steady-state solution develops an instability for the inten- 
sities of the waves exceeding a threshold |23j . The above 
simple examples suggest that the unstable regimes are 
fundamental in nonlinear optics and this appears to be 
indeed true: instabilities, self-oscillations, pattern forma- 
tion, and spatio-temporal chaos are observable in many 



nonlinear optical systems (see, e.g., Refs. [24|]-[2(| for re- 
views). Even incoherent light beams have been recently 
shown to exhibit pattern formation and 'optical turbu- 
lence' due to the so-called modulation instability [27I . 
indicating that unstable behavior does not require the 
perfect coherence of the underlying wave field. 

IV. ANGULAR CORRELATIONS AND 
COHERENT BACKSCATTERING FOR WAVES 
IN NONLINEAR DISORDERED MEDIA 

After having briefly reviewed the key features of wave 
propagation in linear disordered and homogeneous non- 
linear media, we are now in a position to attack the cen- 
tral problem of the present paper: diffuse wave prop- 
agation in nonlinear disordered media. A particularly 
pictorial description of the latter can be developed if the 
nonlinearity is assumed to be weak enough to have no 
significant effect on the general diffuse character of wave 
propagation in the medium and the values of the diffu- 
sion constant D and the mean free path I. This requires 
the 'nonlinear scattering length' ^nl due to the nonlin- 
ear term in Eq. to be much larger than I. £nl can 
be estimated by assuming the short-range intensity cor- 
relation function in a nonlinear disordered medium to 
be approximately the same is in the linear one [see Eq. 
(J3J)]. Calculating the total scattering crossection of sta- 
tistically isotropic fluctuations of the nonlinear dielectric 
constant e 2 8I (r) in the Born approximation, we obtain 

^NL K k / $(K)KdK, (12) 
JO 

where 

$ (K)«4/(«(r)«(r + Ar))e x p(- i KA^A r ,(13) 

and hence £nl ~ (An 2 /co) _1 , where An = 712/0 is the 
typical value of the nonlinear correction to the refractive 
index, n 2 = £2/2, and I is the typical value of (I(r,t)). 
The condition of weak nonlinearity, I <C ^nl, then be- 
comes 

An 2 fc ^ < 1, (14) 

where, we recall, ^> 1. We adopt this condition 
throughout the rest of the paper. 

We start the analysis of the effect of nonlinearity on 
multiple scattering of waves in disordered media by con- 
sidering the angular correlation of scattered wave fields, 
Caba'b', as defined in Eq. (J5J (Sec. [HJ, assuming the two 
incident waves to have different amplitudes (A for the 
wave with the wave vector k a and A' for the wave with 
the wave vector k^). Due to the nonlinearity of the 
medium, C a ba'b' will now depend not only on k a , k^, 
k;,, and k^, but also on A and A'. For simplicity, we 
put A' — ► (so that the effect of the nonlinearity on the 
propagation of the second wave is negligible) and limit 
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ourselves to the case of instantaneous (tnl = 0) local 
( a NL = 0) nonlincarity. 

The physical origin of the nonlinearity-induced modi- 
fication of Caba'b' can be easily understood in the frame- 
work of the path-integral picture of wave propagation. 
The wave of amplitude A, propagating along some wave 
path of length s through a nonlinear disordered medium, 
acquires an additional, 'nonlinear' phase shift A^nl with 
respect to the wave of amplitude A' — > 0, propagating in 
an essentially linear medium: 



A0nl = k n 2 / rfsiJ(ri), 
Jo 

the average value of the latter being 



(A0 NL ) k n 2 A 2 s, 



where we assume (I(ri)) A 2 . In a linear medium, 
a similar phase difference arises between two waves at 
different frequencies u> ^ uj': 



A0aw = (Aw/c)s, 



(17) 



where Aoj = lo — lo' <C lu. It is well-known that the latter 
dephasing modifies the angular correlation function of 
scattered fields, which becomes (in transmission) Q 



a 



ba'b' 



= 5 



Aq„,Aq 



(Aq 2 a - 2z 7 2 ) 1 / 2 il , (18) 



where -y 2 = 3Aw/(2£c). 

The formal analogy between Eqs. l(T^|l and l(T7)l sug- 
gests that as long as the nonlinearity affects only the 
phase of the wave and not its amplitude (i.e. for a weak 
nonlinearity), the correlation function C a ba'b' hi a non- 
linear disordered medium should be approximately given 
by Eq. I|18l) . where Aw is replaced by cun 2 A 2 . This is 
indeed confirmed by a perturbative diagrammatic cal- 
culation 01 [ see t ne diagram of Fig. ^a)], that allows 
one to obtain the lowest term in the series expansion of 
Caba'b' hi n 2 A 2 <C 1. Applying then a 'natural' ansatz 
f(x) - [ie 2 /{2x)]f'(x) ~ f[(x 2 - te 2 ) 1 ' 2 ] yields 



Ca- 



ba'b' 



<W,A qb F T { [Aq 2 a L 2 2l{L/tf] } , (19) 



where £ is a new characteristic length, characterizing the 
damping of correlation due to the nonlinear effects: 



koAn 



(20) 



Here An w n 2 A 2 is the typical value of the nonlinear 
part of the refractive index. The physical meaning of 
£ is clear: it is the characteristic length for the loss of 
phase coherence between two waves following the same 
diffusion path, one with a vanishing amplitude A' — » 
and the other one with a finite amplitude A. By analogy 
with the case of frequency correlation [Eq. ljT%|) ]. 7nl can 
be defined as 7nl = !/£• The short-range correlation 




(15) (b) 



(16) 



c i ;> 




FIG. 1: (a) The diagram contributing to the angular correla- 
tion function of scattered wave fields in a nonlinear disordered 
medium, (b) Short-range intensity correlation function for a 
diffuse wave transmitted through a slab. The figures are from 
Ref. [17]. 



function of intensity fluctuations, , can be obtained 
by taking a square of the absolute value of Eq. (|19|l . We 
show this correlation function in Fig.^b). 

Similarly to the case of transmission geometry, corre- 
lation between diffusely reflected waves can also be cal- 
culated taking into account the additional dephasing due 
to the nonlinearity of the medium. Again, assuming that 
the two incident waves have different amplitudes A and 
A' — > 0, one finds for a semi-infinite disordered medium 



Caba'b' * S AclaAclb {Fr [(Aq 2 a - 2i(l + 0/e) 1 



1/2 



q a + q h + Aq a | 2 -2i(l + C)/e 



1/2' 



(21) 



This result is to be compared to the 'linear' result (7J|. 
In Fig. Efa) we show the experimentally accessible short- 
range correlation function of intensity fluctuations, C^ , 
obtained by taking a square of the absolute value of Eq. 
(|21ll . Just as in a linear medium, two peaks occur at 
Aq a = and q a + q/, + Aq a = 0, respectively. 

Finally, the coherent backscattering cone can be cal- 
culated at the same level of approximation ^tJ- If t ne 
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(a) 




very beginning [Eq. i|16|) ]. The intensity fluctuations in 
disordered media are, however, known to be important 
and exhibit some nontrivial features [such as, e.g., long- 
range correlations, see Eq. Q]. It is therefore important 
to analyze the role of these fluctuations in the presence 
of nonlinearity, when the highly irregular spatial inten- 
sity distribution (speckle pattern) can give rise to a ran- 
dom modification of the dielectric constant (and, conse- 
quently, refractive index). Such an analysis is presented 
in the following section. 



(b) 




INSTABILITY OF DIFFUSE WAVES IN 
NONLINEAR DISORDERED MEDIA 



V. 



As we already mentioned in Sec. IIIII unstable regimes 
are encountered in many nonlinear optical systems. It 
is therefore natural to expect that they may appear in a 
disordered nonlinear medium as well. The origin of the 
instability in this case is relatively easy to understand 
in the framework of the path-integral picture of wave 
propagation [T^ . 



FIG. 2: (a) Short-range angular correlation function of the 
fluctuations of diffusely reflected intensity in a semi-infinite 
nonlinear disordered medium, (b) Line shape of the coher- 
ent backscattering cone in a semi-infinite nonlinear disordered 
medium for L — > co [vanishing linear and nonlinear absorp- 
tion, curve (a)] and I/ a L = W£ [curve (b)]. The figures are 
from Ref. [17]. 



nonlinear coefficient £2 in Eq. is purely real, no devia- 
tion from the linear result is found within the present 
theoretical framework. The contribution of the crossed 
diagrams to the intensity of backscattered wave is, how- 
ever, modified if £2 has a small imaginary part (i.e. in 
the presence of nonlinear absorption). This phenomenon 
can be described by introducing a generalized absorption 
length £^ L , accounting for both linear and nonlinear ab- 
sorption: 

(Z^ ?s Z| + T (1 + c)A2lme2 - (22) 

Additional absorption, introduced by the nonlinear term 
in Eq. will lead to the rounding of the triangular 
peak shape of the coherent backscattering cone, as illus- 
trated in Fig. Efb) . Mathematically, this phenomenon is 
described by Eq. © with L a replaced by L^ h . 

To conclude this section, we would like to emphasize 
that in calculating the angular correlation functions [Eqs. 
iJTDJ, JjHJ and Figs. Gib), Ota)] and the coherent backscat- 
tering cone [Fig- Etb)] , we have completely neglected the 
fluctuations of the wave intensity inside the disordered 
medium, replacing J(r) by (I{r)) ~ A 2 = const in the 



A. Heuristic description 

In the path-integral picture, the diffusely scattered 
wave field at some position r inside the disordered 
medium is represented as a sum of partial waves trav- 
eling along various diffuse paths. Due to the nonlinearity 
of the medium, the phase <p(r, t) of a given partial wave 
arriving at r at time t depends on the spatio-temporal 
distribution of intensity 7(r, t) (speckle pattern) inside 
the medium. At the same time, J(r, t) is a result of in- 
terference of many partial waves, and, therefore, I(r, t) 
is sensitive to their phases </>(r, t) . This leads to a sort of 
feedback mechanism: The phases of partial waves depend 
on the speckle pattern I(r, t), while the latter depends on 
the phases. It appears that such a feedback can desta- 
bilize the time- independent speckle pattern I(r), leading 
to 7(r, t) that fluctuates spontaneously with time, if the 
nonlinearity is strong enough [16j. As we show below, 
this phenomenon can be understood in reasonably sim- 
ple terms in the case of local (onl = 0) and instantaneous 
(tnl = 0) Kerr nonlinearity. 

Consider a partial wave traveling along a typical diffuse 
path of length so(t) ~ L 2 /£ from the source of waves to 
some point r located at a distance L from the source. 
We assume that the path length So(t) can vary slowly 
with time due to, e.g., the motion of scattering centers 
[corresponding to a time-dependent disorder Se(r, t) in 
Eq. JJJ], while the initial and final points of the path are 
fixed. The variations of sq with time are assumed to be 
slow (i.e. so does not change significantly during the time 
Td = L 2 / D required for the wave to cover the distance of 
order sq). The phase <fi(r, t) of the wave contains a 'linear' 
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contribution </>L(r, t) = koSo(t) and the 'nonlinear' one 

r s o{t) 

<fah(r,t) = k dsi[Ae NL (ri,ii)/2] 
Jo 

r s o(t) 

= k n 2 / dsil(n,ti), (23) 
Jo 

where A£nl(i"i,£i)/2 = n 2 I(ri,t\) is the nonlinear part 
of the refractive index, s± is a curvilinear coordinate of 
the point ri along the path, t% = t — [so(t) — si]/c is the 
time of wave passage through ri, and the integrals are 
along the path. 

Let the fluctuating part of the linear dielectric constant 
Se(r, t) to change slightly during some time interval At: 
5e(r,t) -» 5e(r,t + At) = Se(r,t) + A[Se(r,t)], leading 
to corresponding changes A/(r, t, At) and A<fr(r, t, At) 
of intensity and phase, respectively. The latter can be, 
in principle, found from Eq. A(j> contains the linear 
and nonlinear parts, A0l and A0nl, respectively. Under 
very general assumptions, one can show that (AI) = 
and (Atfi) = 0, where the averaging is over 5e(r,t), 
A[fe(r, t)), and over all paths of the same length sq. The 
second moment (A</> 2 ) = (A^) + (A<^ L ) > 0, where 
(A^l) is calculated in the framework of the so-called 
diffusing-wave spectroscopy (DWS) [5j and depends on 
the way in which 5e(r,t) is modified. If, for exam- 
ple, the disordered medium is a suspension of Brown- 
ian particles (particle diffusion coefficient Db), (A^l) = 
(At/T )(s /£), where r = {A^Db)' 1 Q. For the vari- 
ance of the nonlinear phase difference, Eq. (|23[1 yields 

/■so rso 
(A</>nl) = ^0 n 2 / ds l / ds 2 

Jo Jo 
(AJfn.tOA/Cra.ta)), (24) 

where t 2 = f i + At and both integrals are along the same 
diffusion path. Obviously, 

(A/(r 1 ,t 1 )A/(r 2 ,i 2 ))=2[(«(r 1 ,^/(r 2 ,t)) 
- (M(n,t)*J(r 2 ,t + At))] 
~ (AcPt) (51(^,1)51(^,1)} , (25) 

where the second line is obtained by assuming that 
(<5/(ri, t)5I(r 2 , t + At)) in the nonlinear medium is close 
to its value in the linear one (a sort of perturbation the- 
ory), replacing (5I(ri,t)5I(r 2 ,t + At)) by the 'linear' re- 
sult: 

(5I(T 1 ,t)SI(r 2 ,t + At)) ~ (5I(v u t)8I{r 2 ,t)) 

= exp [-(1/2) (A^>] , (26) 

and assuming (A(f>1) <C 1. 

Substituting Eq. lt25|) into Eq. (PU) and changing the 
variables of integration to s = (s\ + s 2 )/2 and As = 
si — s 2 , we obtain 

(A4 L ) = (Atfc) k 2 n 2 2 / ds / dAs (51(^)51^)) . 

Jo J-s 

(27) 



As we discussed in Sec. |H1 the correlation function of 
intensity fluctuations (5I(ri)5I(r 2 )) in Eq. I|27|) contains 
two contributions: a short-range one [Eq. ®] and a long- 
range one [Eq. Q]. To perform the integration in Eq. 
(|27|l . we replace Ar = \y\ — r 2 \ in the expression (J3J) 
for the short-range correlation function by As, assuming 
the wave path to be ballistic at distances shorter than I. 
In contrast, for Ar > £ the wave path is diffusive, and 
hence Ar in the expression Q for the long-range correla- 
tion function can be substituted by {Asi) 1 / 2 . Performing 
then integrations in Eq. (|2"7J) yields 

(A^ L >~p(A^), (28) 
where we introduce the bifurcation parameter [T^ | 

p = An 2 (j\ (k t+j\ (29) 

and the numerical factors of order unity are omitted. 

Note that Eq. (|29|l is a sum of two contributions. The 
first one, An 2 (L/£) 2 ko£, originates from the short-range 
correlation of intensity fluctuations [Eq. while the 
second contribution, An 2 (L/£) 3 , is due to the long-range 
correlation [Eq. If p <C 1, the nonlinear term 

(A0^ L ) represents just a small correction to the lin- 
ear one (A0l) and our perturbation approach is likely 
to be valid. In contrast, the above perturbation theory 
diverges for p > 1, since (A^l) becomes larger than 
(A^l)- This suggests that p ~ 1 is a critical point be- 
yond which (i.e. for p > 1) the physics of the nonlin- 
ear problem is no longer similar to the physics of the 
linear one and hence the perturbation theory cannot be 
used. Although the mere breakdown of the perturbation 
theory does not permit to draw any far-reaching conclu- 
sions about the speckle pattern beyond the critical point 
p ~ 1, a more rigorous analysis (see Sees. IV 51 and IVT^) 
shows that p ~ 1 defines the instability threshold of the 
multiple-scattering speckle pattern, and that at p > 1 
the latter should exhibit spontaneous temporal fluctua- 
tions even for a time-independent disorder 5s(r) (i.e. for 
vanishing mobility of scattering centers: Db —> and 
To — > oo in the case of Brownian motion) . 

It is pertinent to note the extensive nature of the bi- 
furcation parameter p. According to Eq. I|29l) , weak non- 
linearity can be efficiently compensated by a sufficiently 
large sample size L and p ~ 1 can be reached even at van- 
ishing An, provided that L is large enough. In the limit 
of large sample size L ^> ko£ 2 , one gets p ~ An 2 (L/£) 3 
(as also found in Ref. 28] ). This result is due to the 
long-range correlation Q of intensity fluctuations which 
dominates in the integral of Eq. 1)2 7|l because of its slow 
decrease (cx |ri — r 2 \ ) and consequent extension over 
the whole disordered sample. Another remarkable fea- 
ture of Eq. H29fl is the independence of p of the sign of 
An. The phenomenon of the speckle pattern instabil- 
ity is therefore expected to develop in a similar way for 
both positive and negative nonlinear coefficients n 2 . This 
is not common for the instabilities in nonlinear systems 
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without disorder [^-[^l since the latter are often re- 
lated to the self-focusing effect, arising at n 2 > only. 
The instability discussed in the present paper is of dif- 
ferent type and has nothing to do with the self-focusing 
[which is negligible if the condition l|14l) is satisfied] . This 
does not mean, however, that the phenomena similar to 
those discussed in the present paper do not exist in ho- 
mogeneous nonlinear media. In fact, it is easy to see 
that E g. (12311 describes nothing but the self-phase modu- 
lation [13 of the multiple-scattering speckle pattern. It 
is worth noting that the development of the self-phase 
modulation in a disordered medium appears to be rather 
different from that in the homogeneous case. Indeed, in a 
homogeneous medium of size L the nonlinear phase shift 
is deterministic and </>nl = kon2loL for a wave of inten- 
sity Iq, while in the case of diffuse multiple scattering </>nl 
is random with the average value (0nl) — kon2lo(L 2 /£) 
and variance (</>nl) — (0nl) — P- The threshold of the 
speckle pattern instability p ~ 1 is then simply the point 
where the sample-to-sample fluctuations of the nonlinear 
phase shift become of order unity. 

To conclude this subsection, we would like to comment 
on the recent results on modulation instability (MI) of 
incoherent light beams in homogeneous nonlinear me- 
dia [23. In this 2D speckle pattern (speckled 
beam) with a controlled degree of spatial and temporal 
coherence is 'prepared' by sending a coherent laser beam 
on a rotating diffuser. The beam is then incident on a 
homogeneous nonlinear medium (inorganic photorefrac- 
tive crystal). Pattern formation and 'optical turbulence' 
is observed if the intensity of the beam exceeds a spe- 
cific threshold, determined by the degree of coherence 
of the beam (coherent MI has no threshold). Although 
the study presented in this paper is also concerned with 
the instability of speckle patterns in nonlinear media, it 
differs essentially from that of Ref. [2]j at least in the fol- 
lowing important aspects: (a) the initially coherent wave 
loses its spatial coherence inside the medium, due to the 
multiple scattering on heterogeneities of the refractive in- 
dex, (b) the temporal coherence of the incident light is 
perfect (in the case of incoherent MI [27J, the coherence 
time of the incident beam is much shorter than the re- 
sponse time of the medium tnl), and (c) the incident 
light beam is completely destroyed by the scattering af- 
ter a distance ~ £ and light propagation is diffusive in the 
bulk of the sample. In addition, diffuse multiple scatter- 
ing of light results in a distributed feedback mechanism, 
absent in the case of MI. 



B. Path-Integral Picture 

A mathematically rigorous formulation of the heuristic 
analysis presented in the previous subsection has been 
performed in Ref. 0]. In addition to the condition of 
weak nonlinearity Ijl4|l we anticipate that the expected 
spontaneous fluctuations of the speckle pattern beyond 
the instability threshold will be slow enough or, more 
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FIG. 3: (a) Graphical solution of Eq. (J3UJ at k £ = 100, 
I An| = 10 -4 and the values of £/ indicated near each curve. 
The solid lines show F(/3), the dashed line is exp( — f3). Thick 
solid lines with arrows illustrate the iterative solution of Eq. 
13011 . (b) Bifurcation diagram of the speckle pattern in a 
semi-infinite disordered medium for ko£ = 100 and the values 
of 1/ indicated near the each curve. The threshold values 
of |An|, following from the condition p = 1 at £/L R > 0, are 
indicated by arrows. In the absence of absorption, £/L a — 
(dashed line), there is no threshold and the speckle pattern is 
unstable at any (even infinitely small !) [An|. The figures are 
from Ref. [16]. 



precisely, that the typical time of the speckle pattern dy- 
namics will be much longer than the typical time between 
two successive scattering events £/c. 

For the sake of the mathematical simplicity, we con- 
sider a plane monochromatic wave incident on the sta- 
tionary [fe(r) is time- independent] semi-infinite disor- 
dered medium (L —* 00) with a finite (macroscopic) 
absorption length L a ^ £ and a local (onl = 0) 
instantaneous (tnl = 0) Kerr nonlinearity. It ap- 
pears that instead of replacing the intensity correlation 
function (<5J(ri, t)SI(r 2 , t + At)) in Eq. (J23) by its 'lin- 
ear' value (and thus constructing a sort of perturba- 
tion theory), one can obtain a self-consistent descrip- 
tion of the speckle instability. The analysis [l6j leads 
to a self-consistent equation for [3 = — In g\ , where 

9l ( T ) = (E(t)E*(t + T)} / (^\E(t)\ 2 ^ is the time autocor- 
relation function of diffusely reflected wave, and g\ de- 
notes lim-r^oo gi (t): 

exp(-/?)=F(/3), (30) 

where F{(3) is a known monotonically decaying function 
and F(0) = 1. 

Fig. OJa) illustrates the graphical solution of Eq. 1|3U|) 
for some realistic values of parameters. At strong enough 
absorption [two upper curves in Fig.^a)], a unique solu- 
tion (3 = Pi = exists, corresponding to g\ — exp(— (3) = 
1. This means that the time autocorrelation function 
of diffusely reflected wave g\ (r) does not decay with 
time, remaining equal to 1 for all t (even for t — * 00). 
The speckle pattern is therefore stationary (i.e. docs 
not change with time). In contrast, for weak absorp- 
tion [three lower curves in Fig. GJa)] a second solution 
(3 = P2 > arises, corresponding to gi = exp(— [3) < 1. 
Eq. i|3Ufl now has two fixed points: /3 = j3\ = and 
j3 = /fe > 0. It is seen from Fig. EUa) (see thick solid 
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lines with arrows) that an iterative solution of Eq. i|30|) 
starting at some (3 > converges to (3 = f3\ = for 
strong enough absorption and to (3 = 02 > for weak 
absorption. (3 = f32 > is therefore a stationary point of 
Eq. I|30l) in the latter case, corresponding to the physi- 
cally realizable solution beyond the instability threshold. 
(3 > and g\ < 1 signifies that the time autocorrela- 
tion function of diffusely reflected wave g± (r) decays with 
r despite our assumption of stationary disorder [time- 
independent Se(r)], although the present analysis does 
not allow us to estimate the characteristic time scale of 
this decrease. Decaying <?i(r) corresponds to an unstable 
(i.e. time-varying) speckle pattern. Note that in con- 
trast to the fluctuations of the speckle pattern due to the 
motion of scattering centers, employed in DWS to study 
the dynamics of the medium |5| , the speckle dynamics in 
a stationary nonlinear random medium is spontaneous, 
i.e. it does not originate from the motion of scattering 
centers (since the latter are immobile) but is intrinsic for 
the underlying nonlinear wave equation . Spontaneous 
speckle dynamics is irreversible and hence the instability 
threshold p ~ 1 is the point where the time-reversal sym- 
metry is spontaneously broken for the multiple-scattered 
waves. 

Fig- Bib) shows the time autocorrelation function g\ — 
cxp(— (3) obtained by solving Eq. (|30|l at k$£ — 100 and 
several values of £/L & . The transition from the stable 
(gi = 1) to unstable (<7i < 1) regime is clearly seen when 
£/L a > 0. At £/L a = (no absorption), we find that 
g% < 1 at any, even infinitely small |An|. The border 
between stable (g\ = 1) and unstable (gi < 1) regimes 
can be found analytically by requiring (d/d(3)F((3) = —1 
at (3 = 0. This yields p ~ 1 as the instability threshold, 
where the bifurcation parameter p is defined by Eq. i|29|) 
with L replaced by L a . The threshold values of |An| 
following from the condition p = 1 are shown in Fig. 
Bib) by arrows. 



C. Langevin Description 

Although the path-integral approach, sketched in the 
two previous subsections, allows us to predict the insta- 
bility of the multiple-scattering speckle pattern in a non- 
linear disordered medium and even to calculate g\{r) in 
the limit r — > oo, it does not permit to estimate the time 
scale of spontaneous intensity fluctuations beyond the in- 
stability threshold. Meanwhile, this issue is of primary 
importance in view of the possible experimental observa- 
tion of the instability phenomenon. In this subsection, 
we show that the spontaneous speckle d yna mics can be 
studied by using the Langevin approach [I^, [ltj • 

As in the previous subsections, we assume the non- 
linearity to be local (onl = 0), but consider arbitrary 
nonlinearity response time tnl- The nonlinear part of 
the dielectric constant is assumed to be governed by Eq. 
IjlUp. In addition, we neglect the absorption (assuming 
L a 3> L) and restrict our consideration to the case when 



the development of the speckle pattern instability is dom- 
inated by the long-range intensity correlations. This re- 
quires the sample size to be large enough [L » ko£ 2 , 
see Eq. (|29[1 and its accompanying discussion] and the 
speckle dynamics to be slow (time scale of spontaneous 
intensity fluctuations r > T D [(k £ 2 )/ L] 2 ). The latter 
condition follows from Eq. i|27[l . where the limits of in- 
tegration over As should be set to =Fct, if cr <C sq (or, 
equivalently, r <C To), since the correlation has simply 
no time to establish beyond these limits due to the finite 
speed of wave propagation. Integration in Eq. l|27l) then 
yields 



(Acf 



2 \ 

NL/ 



An' 



k l- 



L 
1 



1/2' 



(31) 



where the first term in square brackets [the same as in Eq. 
(|29|l ] results from the short-range correlation J2J), while 
the second term is due to the long-range one l@J ■ Requir- 
ing that the second term dominates the first, we obtain 
t > T u [(k Q £ 2 ) / L} 2 as stated above. 

The Langevin equation for the intensity fluctuation 
5I(r,t) reads 0: 



£V 2 <5/(r,i) = -V-j oxt (r,i), 



(32) 



where j e xt(r, t) are random external Langevin cur- 
rents that have zero mean and a correlation function 
(jixt(M)ii x t( r i>*l)) given by the diagram (i) of Fig. 
01a). Being a 'fingerprint' of disorder, the Langevin cur- 
rents will be modified if we modify the dielectric con- 
stant of the medium. In our case, the linear part of 
the dielectric constant is fixed [Se(r) in Eq. is time- 
independent], but its nonlinear part Ae^^^t) can vary 
with time, and it is relatively easy to obtain the following 
dynamic equation for j cx t(r, t) 



isy to oc 



d_ 

di 



jcxt(r,i) 



d 3 r' / dAt 
v Jo 



q(r,r',At) — Ae NL (r',i - At), (33) 

where q(r,r',A£) is a random response func- 
tion with zero mean and the correlation function 
(g^( r : r '> At)q ( -^*(r 1 , r[, Ai x )) given by a sum of the 
diagrams (ii) and (iii) of Fig. Ela). Eqs. (f32J an( l EH 
together with Eq. I|1U|) for Aenl form a self-consistent 
system of equations for the stability analysis of the 
speckle pattern. 

Consider now an infinitesimal periodic excitation of the 
stationary speckle pattern: SI(r,t) — SI(r,a) exp(ai), 
where a = iQ + A ^ and > 0. Such an excitation can 
be either damped or amplified, depending on the sign of 
the Lyapunov exponent A. The value of A is determined 
by two competing processes: on the one hand, diffusion 
tends to smear the excitation out. while on the other 
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FIG. 4: (a) Diagrams contributing to the correlation func- 
tions of Langevin currents j e xt(r, t) [diagram (i)] and random 
response functions q(r, r',Ai) [diagrams (ii) and (in)]. The 
diagrams (ii) and (iii) are obtained by the functional differen- 
tiation of the diagram (i) with respect to the dielectric con- 
stant of the disordered medium. Wavy lines in the diagrams 
(ii) and (iii) denote fco factors, (b) Surface describing the 
stability of the multiple-scattering speckle pattern in a disor- 
dered medium with fast nonlinearity (Odtnl <C 1). For given 
frequency Q and bifurcation parameter p, the surface shown 
in the figure allows one to determine the value of the Lya- 
punov exponent A. If A > 0, the speckle pattern is unstable 
with respect to excitations at frequency Q. The figures are 
from Ref. [18]. 



hand, the distributed feedback sustains its existence. The 
mathematical description of this competition is provided 
by Eqs. J22), and {TH)) that after the substitution 

of SI(r, t) — 6I(r, a) exp(iai) [and similarly for j ex t(r, t) 
and Aent/f, t)] lead to the following equation for p, fl, 
and A 



■£NT/r,f) 

mm- 

p~F l (ft/Q D ,A/Q D )F 2 (Ot nl ,At nl ) , 



(34) 



where VL-q = 1/Td = D/L 2 and the function F\ is shown 
in Fig. Efb) for the case of the sample with open bound- 
aries, while F 2 (x, y) = 1 + x 2 + y 2 + 2y. 

Let us first consider the fast nonlinearity, assum- 
ing tnl^d <C 1 |lq . In this case, the nonlinear re- 
sponse takes much less time than the typical perturba- 
tion 5I(r, a) needs to extend throughout the disordered 
sample and we can set F 2 (Otnl, Atnl) — 1 in Eq. ijSljl. 
The stability of the speckle pattern with respect to pe- 
riodic excitations is then described by Fig. BJb). For a 
given frequency f2, the sign of the Lyapunov exponent A 
depends on the value of p. Excitations at frequencies O 
corresponding to A < are damped exponentially and 
thus soon disappear. In contrast, excitations at frequen- 
cies Q corresponding to A > are exponentially ampli- 
fied, which signifies the instability of the speckle pattern 
with respect to excitations at such frequencies. Noting 
that A is always negative for p < 1, we conclude that 
all excitation are damped in this case and the speckle 
pattern is absolutely stable. In an experiment, any spon- 
taneous excitation of the static speckle pattern will be 
suppressed and the speckle pattern will be independent of 
time: <5/(r, t) — 51 (r), as in the linear case. When p > 1, 
an interval of frequencies < Q < O max starts to open 
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FIG. 5: (a) Frequency-dependent instability threshold for a 
disordered sample with instantaneous nonlinear response and 
open boundaries. The inset shows the maximal Lyapunov 
exponent as a function of the bifurcation parameter p. (b) 
The same as (a) but for a sample with a nonzero nonlinearity 
response time tnl. The dashed line shows the result corre- 
sponding to tnl^d = [the same as the solid line in the panel 
(a)]. The figures are from Ref. [18] (a) and [19] (b). 



up with A > 0. The speckle pattern thus becomes un- 
stable with respect to excitations at low frequencies. In 
an experiment, any spontaneous excitation of the static 
speckle pattern at frequency Q € (0, f2 max ) will be am- 
plified and one will observe a time-varying speckle pat- 
tern 5I(r, t). Note that the absolute instability threshold 
p ~ 1 agrees with the results of Sees. IV Al and IV hfl while 
completely different calculation techniques have been ap- 
plied in the three cases. 

Fig. UJa) shows the frequency-dependent instability 
threshold for the disordered sample with open bound- 
aries and fast nonlinear response. Detailed analysis 0] 
of Eq. H34f) shows that the frequency-dependent thresh- 
old value of p scales as 1 -I- Ci(fl/flv) 2 (where C\ ~ 1 is 
a numerical constant) at fl <C Od and as {Q./Q.-a) 1 / 2 at 
£1 3> £1d • The latter result can also be obtained from Eq. 
Ipnjl by setting r ~ I/O and applying (A0 NL ) > (A0£) 
as the instability condition. For a given value of p > 1, 
the maximum excited frequency is f2 max ~ ^d(p — l) 1 ^ 2 
for p — 1 <C 1 and O max ~ QdP 2 for p — 1 > 1 . 

Fig-EIb) illustrates the effect of noninstantaneous non- 
linearity on the frequency-dependent instability thresh- 
old for several values of tnl^d uM- Obviously, slow 
nonlinear response of the medium rises the instability 
threshold for high-frequency excitations, while the abso- 
lute instability threshold p ~ 1 remains the same, in- 
dependent of the nonlinearity response time tnl- The 
reason for this is that just above p — 1 the speckle pat- 
tern becomes unstable with respect to excitations at very 
low frequencies, for which £1 <C 1/tnl ^ s always fulfilled 
and which, therefore, are not sensitive to the value of 
tnl [mathematically, ^(Htnl, Atnl) — 1 in Eq. 
In the limit of slow nonlinearity (tnl^d S> 1) and for 
1 < p — 1 < 1 we can set F\ (tl/Clo, A/0 D ) ~ 1 in 
Eq. (|34|) and find the threshold value of p to scale as 
1 + (Ot nl ) 2 and il max ~ T^l(p - I) 1 ' 2 . 

As follows from the above analysis, a continuous low- 
frequency spectrum of frequencies (0, O max ) is excited at 
p > 1. In addition, the Lyapunov exponent A appears to 
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decrease monotonically with £1, favoring no specific fre- 
quency fl just above the threshold |l9|. This allows us to 
hypothesize that at p = 1 the speckle pattern undergoes 
a transition from a stationary to chaotic state. Such a 
behavior should be contrasted from the 'route to chaos' 
through a sequence of bifurcations, characteristic of many 
nonlinear and, in particular, optical systems |22^-|2fi|. 

VI. CONCLUSION 

The present paper reviews some of the recent theoret- 
ical developments in the field of multiple wave scattering 
in nonlinear disordered media. To be specific, we consider 
optical waves and restrict ourselves to the case of Kerr 
nonlinearity. Assuming that the nonlinearity is weak, we 
derive the expressions for the angular correlation func- 
tions and the coherent backscattering cone in a nonlinear 
disordered medium fSec. ITvjl . In both transmission and 
reflection, the short-range angular correlation functions 
of intensity fluctuations for two waves with different am- 
plitudes (A and A' — > 0) appear to be given by the same 
expressions [Eqs. © and 0, respectively] as the angu- 
lar correlation functions for waves at two different fre- 
quencies (w and lj' = uj — Au>) in a linear medium, with 
Ato replaced by 2^c/(3£ 2 ), where £ is a new nonlinear 
characteristic length defined by Eq. (|2C)|> . The coherent 
backscattering cone is not affected by the nonlinearity, 
as long as the nonlinear coefficient £2 in Eq. Q is purely 
real. If £2 has an imaginary part (which corresponds 
to the nonlinear absorption), the line shape of the cone 
is given by the same expression @ as in an absorbing 
linear medium, where the linear macroscopic absorption 
length L a should be replaced by the generalized absorp- 
tion length L^ L defined by Eq. 

For the nonlinearity strength exceeding a threshold 
p ~ 1 [with the bifurcation parameter p given by Eq. 
(|29[) ]. we predict a new phenomenon — temporal insta- 
bility of the multiple-scattering speckle pattern — to take 
place (Sec. 0). The instability is due to a combined ef- 
fect of the nonlinear self-phase modulation and the dis- 
tributed feedback mechanism provided by multiple scat- 
tering and should manifest itself in spontaneous fluctua- 
tions of the speckle pattern with time. Since the sponta- 
neous dynamics of the speckle pattern is irreversible, the 
time-reversal symmetry is spontaneously broken when p 
surpasses 1. The important feature of our result is the 
extensive nature of the instability threshold, leading to 
an interesting possibility of obtaining unstable regimes 



even at very weak nonlinearities, provided that the dis- 
ordered sample is large enough. To study the dynamics of 
multiple-scattering speckle patterns beyond the instabil- 
ity threshold, we generalize the Langevin description of 
wave diffusion in disordered media fSec. IV Explicit 
expressions for the characteristic time scale of sponta- 
neous intensity fluctuations are derived with account for 
the noninstantaneous nature of the nonlinearity. The 
results of this study allow us to hypothesize that the dy- 
namics of the speckle pattern may become chaotic im- 
mediately beyond the instability threshold, and that the 
cascade of bifurcations, typical for chaotic transitions in 
many known nonlinear systems |22|-|26||. might not be 
present in the considered case of diffuse waves. 

Finally, we discuss the experimental implications of our 
results. With common nonlinear materials, such as, e.g., 
carbon disulfide, ri2 ~ 3 x 10~ 14 cm 2 /W can be realized 
P3- At / ~ 1 MW/cm 2 this yields An ~ 3 x 10" 8 and 
the characteristic length £ defined in Eq. 120(1 is £ ~ 2 cm 
for £ ~ 100 /jm and A ~ 0.5 /im. As follows from Eq. lTH?)l . 
£ < L is required to observe a sizeable effect of the nonlin- 
earity on the angular correlation function of transmitted 
light, and hence using a 2 cm-thick disordered sample 
should suffice to make the effect of nonlinearity measur- 
able in, e.g., a dense suspension of carbon disulfide parti- 
cles. In contrast, £ < £ is necessary to observe the effect 
of the nonlinearity on the angular correlation function in 
the reflection geometry. This requires a stronger nonlin- 
earity. In nematic liquid crystals, for example, An up to 
0.1 is achievable (see, e.g., Ref. Taking An = 10~ 3 

and I ~ 1 mm [3J|, we get £ ~ 1 mm ~ I. In the case of 
the coherent backscattering cone, very accurate experi- 
mental techniques developed recently to measure the an- 
gular dependence of backscattered light give a hope 
that the effect of the nonlinearity can be observable with- 
out any particular problem. Finally, observation of the 
temporal instability of multiple-scattering speckle pat- 
terns, predicted in Sec. will require, first of all, large 
enough sample size L and as low absorption as possi- 
ble. In the absence of absorption (or for the macroscopic 
absorption length L a > L \ AnNL ~ 10~ 2 (realistic in 
nematic liquid crystals 29]) and L/£ ~ 20 will suffice to 
get p ~ 1 and reach the instability threshold. We note 
that in liquid crystals, the nonlinearity is due to the reori- 
entation of molecules under the influence of the electric 
field of the electromagnetic wave and hence is essentially 
noninstantaneous. This emphasizes the importance of 
including the noninstantaneous nature of nonlinear re- 
sponse in our analysis fSec, IVC|l . 
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